Accelerating a 3D finite-difference wave propagation code using GPU graphics cards

International audienceWe accelerate a three-dimensional finite-difference in the time domain (FDTD) wave propagation code by a factor between about 20 and 60 compared to a serial implementation using Graphics Processing Unit (GPU) computing on NVIDIA graphics cards with the CUDA programming language. We describe the implementation of the code in CUDA to simulate the propagation of seismic waves in a heterogeneous elastic medium. We also implement Convolution Perfectly Matched Layers (CPMLs) on the graphics cards to efficiently absorb outgoing waves on the fictitious edges of the grid. We show that the code that runs on a graphics card gives the expected results by comparing our results to those obtained by running the same simulation on a classical processor core. The methodology that we present can be used for Maxwell's equations as well because their form is similar to that of the seismic wave equation written in velocity vector and stress tensor

Collaborative simulation and scientific big data analysis: Illustration for sustainability in natural hazards management and chemical process engineering

Classical approaches for remote visualization and collaboration used in Computer-Aided Design and Engineering (CAD/E) applications are no longer appropriate due to the increasing amount of data generated, especially using standard networks. We introduce a lightweight and computing platform for scientific simulation, collaboration in engineering, 3D visualization and big data management. This ICT based platform provides scientists an “easy-to-integrate” generic tool, thus enabling worldwide collaboration and remote processing for any kind of data. The service-oriented architecture is based on the cloud computing paradigm and relies on standard internet technologies to be efficient on a large panel of networks and clients. In this paper, we discuss the need of innovations in (i) pre and post processing visualization services, (ii) 3D large scientific data set scalable compression and transmission methods, (iii) collaborative virtual environments, and (iv) collaboration in multi-domains of CAD/E. We propose our open platform for collaborative simulation and scientific big data analysis. This platform is now available as an open project with all core components licensed under LGPL V2.1. We provide two examples of usage of the platform in CAD/E for sustainability engineering from one academic application and one industrial case study. Firstly, we consider chemical process engineering showing the development of a domain specific service. With the rise of global warming issues and with growing importance granted to sustainable development, chemical process engineering has turned to think more and more environmentally. Indeed, the chemical engineer has now taken into account not only the engineering and economic criteria of the process, but also its environmental and social performances. Secondly, an example of natural hazards management illustrates the efficiency of our approach for remote collaboration that involves big data exchange and analysis between distant locations. Finally we underline the platform benefits and we open our platform through next activities in innovation techniques and inventive design

High-frequency simulations of global seismic wave propagation using SPECFEM3D_GLOBE on 62K processors

SPECFEM3D_GLOBE is a spectral element application enabling the simulation of global seismic wave propagation in 3D anelastic, anisotropic, rotating and self-gravitating Earth models at unprecedented resolution. A fundamental challenge in global seismology is to model the propagation of waves with periods between 1 and 2 seconds, the highest frequency signals that can propagate clear across the Earth. These waves help reveal the 3D structure of the Earth's deep interior and can be compared to seismographic recordings. We broke the 2 second barrier using the 62K processor Ranger system at TACC. Indeed we broke the barrier using just half of Ranger, by reaching a period of 1.84 seconds with sustained 28.7 Tflops on 32K processors. We obtained similar results on the XT4 Franklin system at NERSC and the XT4 Kraken system at University of Tennessee Knoxville, while a similar run on the 28K processor Jaguar system at ORNL, which has better memory bandwidth per processor, sustained 35.7 Tflops (a higher flops rate) with a 1.94 shortest period. Thus we have enabled a powerful new tool for seismic wave simulation, one that operates in the same frequency regimes as nature; in seismology there is no need to pursue periods much smaller because higher frequency signals do not propagate across the entire globe. We employed performance modeling methods to identify performance bottlenecks and worked through issues of parallel I/O and scalability. Improved mesh design and numbering results in excellent load balancing and few cache misses. The primary achievements are not just the scalability and high teraflops number, but a historic step towards understanding the physics and chemistry of the Earth's interior at unprecedented resolution

Towards seismic wave modeling on heterogeneous many-core architectures using task-based runtime system

International audienceUnderstanding three-dimensional seismic wave propagation in complex media remains one of the main challenges of quantitative seismology. Because of its simplicity and numerical efficiency, the finite-differences method is one of the standard techniques implemented to consider the elastodynamics equation. Additionally, this class of modelling heavily relies on parallel architectures in order to tackle large scale geometries including a detailed description of the physics. Last decade, significant efforts have been devoted towards efficient implementation of the finite-differences methods on emerging architectures. These contributions have demonstrated their efficiency leading to robust industrial applications. The growing representation of heterogeneous architectures combining general purpose multicore platforms and accelerators leads to redesign current parallel application. In this paper, we consider StarPU task-based runtime system in order to harness the power of heterogeneous CPU+GPU computing nodes. We detail our implementation and compare the performance obtained with the classical CPU or GPU only versions. Preliminary results demonstrate significant speedups in comparison with the best implementation suitable for homogeneous cores

Porting a high-order finite-element earthquake modeling application to NVIDIA graphics cards using CUDA

International audienc

Improved Co-Registration of Sentinel-2 and Landsat-8 Imagery for Earth Surface Motion Measurements

International audienceThe constellation of Landsat-8 and Sentinel-2 optical satellites offers opportunities for a wide range of Earth Observation (EO) applications and scientific studies in Earth sciences mainly related to geohazards. The multi-temporal co-registration accuracy of images provided by both missions is, however, currently not fully satisfactory for change detection, time-series analysis and in particular Earth surface motion measurements. The objective of this work is the development, implementation and test of an automatic processing chain for correcting co-registration artefacts targeting accurate alignment of Sentinel-2 and Landsat-8 imagery for time series analysis. The method relies on dense sub-pixel offset measurements and robust statistics to correct for systematic offsets and striping artefacts. Experimental evaluation at sites with diverse environmental settings is conducted to evaluate the efficiency of the processing chain in comparison with previously proposed routines. The experimental evaluation suggests lower residual offsets than existing methods ranging between RMSE xy = 2.30 and 2.91 m remaining stable for longer time series. A first case study demonstrates the utility of the processor for the monitoring of continuously active landslides. A second case study demonstrates the use of the processor for measuring co-seismic surface displacements indicating an accuracy of 1/5 th of a pixel after corrections and 1/10th of a pixel after calibration with ground measurements. The implemented processing chain is available as an open source tool to support a better exploitation of the growing archives of Sentinel-2 and Landsat-8

High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster

International audienc

Modeling the propagation of elastic waves using spectral elements on a cluster of 192 GPUs

International audienc

High-order Finite-Element Seismic Wave Propagation Modeling with MPI on a large GPU Cluster

International audienc

Accelerating a 3D finite-difference wave propagation code and a spectral-element code using a cluster of GPU graphics cards

International audienc